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Electro-Optical Tactical Decision Aids (TDAs) have proven their utility as tools 
for range performance modeling and mission planning. However, several TDAs are in 
current use in the United States armed forces. In fact, the services use different TDA codes 
which differ in the input data files and their sources required, in the operator expertise 
required, and the hardware required to run the program. 

Within the concept of Joint Operations, which has become crucial in the modern 
battlefield environment, all the services must share procedures, techniques, and often the 
same technology. This thesis presents a comparison between the Army FLER TDA, 
(ACQUIRE), and the infrared module of the Navy/Air Force TDA, WinEOTDA. 
Differences in the modeling of underlying physical principles, input parameters, and 
predicted target detection ranges are presented. Despite differences in input and treatment 
of environmental effects this analysis indicates similar levels of accuracy for the two 
codes. For two scenarios selected average predictions for three “typical” sensors fall 
within 20% of published observations. With further analysis and an operational evaluation 
it may be possible to select one Electro-Optical Tactical Decision Aid for all branches of 
the military. 
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I. INTRODUCTION 



A tactical decision aid (TDA) is defined by Integrated Performance Decisions, 
Inc. at its web page as a tool or set of tools which aid an operator in performing a task. 
Included in this definition are maneuvering recommendations, weapon or decoy 
placement recommendations, target motion analysis, and performance predictions. It 
should help a decision maker (Commander or whoever in charge) by assimilation and 
convenient presentation of data (on targets, own assets and the environment) which of 
themselves are useful, and by providing analysis of the tactical problem beyond what is 
feasible by humans in a timely fashion. 

TDAs can have the form of various types of design product, including 
nomographs, manuals, guides, and of course, TDA programs (codes) resident on desktop 
computers. This thesis focuses on TDA codes for Electro-Optical applications; from now 
on these will be addressed as EOTDA. 

An EOTDA typically consists of input and output processors and three main 
models: a Target Model that determines the inherent signal that emanates from scene 
objects; an Atmospheric Model which estimates the degradation of signal from target to 
sensor due to weather conditions; and a Sensor Model which determines the range at 
which an object can be detected or at which a weapon system can lock-on. 
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In this work, we will evaluate the TDA codes used by the US Navy (EOTDA) and 
US Army (ACQUIRE) for determining and predicting detection ranges, and the 
possibility to use them as standard TDAs for all services or for Joint service use. 

The Army code ACQUIRE is an analytical model based on the static performance 
model FLIR92. It predicts target acquisition performance for systems in which the image 
is in the visible, near infrared, and infrared spectral bands. The detection ranges and 
probabilities predicted by the model represent the expected performance of a set of 
observers (well trained) with respect to an average target having a specified signature and 
size. In this thesis we will discuss the ACQUIRE version 1 model, dated May 1995. 

The Navy code EOTDA is a PC based program, designed mainly as a mission- 
planning tool. It can be used by operational units to predict performance of electro-optical 
weapon systems in the infrared, visible, and laser bands. Prediction is based upon 
meteorological data (forecast), target characteristics, weapon system, and tactics. In this 
thesis we will discuss the EOTDA version 3, dated Jan 1993 and/or the open distribution 
version Win-EOTDA 0.7, dated 1997. 

In this age of smart weapons, more sophisticated methods of weapon delivery 
have been developed. These methods always deal with important and natural phenomena, 
the changing environmental conditions. So, it is necessary to have accurate 
meteorological analyses and forecasts to properly plan and effectively execute tactical 
operations. On the other hand, it is also necessary to present this kind of data in a 
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tactically relevant form in order for it to be easily assimilated for operational purposes. 
Here is where TDAs are very useful as mission-planning tools. 

The first part of this work will deal with a brief discussion of the theory behind 
these codes. Infrared Radiation theory. The next step will be an analysis and comparison 
of how the theory' is implemented in the codes. Finally, a comparison in performance will 
be made of predicted target detection range output from the two programs for equivalent 



operator inputs. 
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II. THEORY AND BACKGROUND 



A. GENERALITIES 



An Electro-optical Tactical Decision Aid (TDAJ is a software model that predicts 
the performance of an air-to-ground weapon system based on environmental and tactical 
information. The two models, ACQUIRE and Win-EOTDA, can perform predictions in 
the visible, near infrared and infrared spectral bands. 

The part of the spectrum that includes our area of interest is shown in Figure 2.1 
that represents the portion of the Infrared (IR) spectrum characterized by atmospheric 
windows; that is, regions in which the atmosphere is transparent to radiation. 

[Ref. 1 :p. 1 79] . 
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Figure 2.1 The Electromagnetic and IR Spectra 



Infrared imaging is the remote sensing and display of infrared flux variations. The 
variations in the displayed image intensity represent apparent temperature variations 
across the scene [Ref.2: p. 7]. A single scene can be represented by an isolated target and 
its surroundings (background). The radiation received from a target is compared with that 
of an equivalent drea of the background. The radiation emitted and reflected from both 
target and background traverse the atmosphere and the aerosol particles suspended in it. 
The target detection is determined by the distribution of temperature over its surface in 
contrast with the temperature distribution of the background, allowing for the attenuating 
effect of the atmospheric path. The next section summarizes information related to 
thermal radiation theory which can be found (and more explicitly detailed) in many 
tutorial sources, for example "The Infrared Handbook" [Ref.3]. 



B. THERMAL RADIATION THEORY 

Optical radiation covers a specific part of the electromagnetic spectrum as we have 
seen before. There are no sharp boundaries between regions; rather, they are distinguished 
from each other by methods of production and detection, transmission through the 
atmosphere, response of the human eye, etc. [Ref. 7:p. 1-25]. 

We start a brief summary of Black Body theory by saying that all objects with 
temperature above absolute zero emit some kind of radiation. A black body is defined as 
an ideal body or surface that absorbs all radiant energy falling upon it, with no reflection 
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or retransmission. As we will see, blackbodies have been studied by many scientists; they 
have stated two important emission properties of a blackbody: 

a. - A blackbody is a perfect emitter of radiation at all angles and all wavelengths. 

b. - The total emitted energy of a blackbody is a function only of its temperature 

[Ref.4:p.45]. 

1. Planck’s Blackbody Law 

Planck hypothesized that a blackbody at a given temperature emits radiation of all 
wavelengths. In particular, the radiant exitance at the various wavelengths integrated over 
a one-micron (jim) band around v, is given by Planck’s law: 

2niiv 5 

Mbb ( v ,t> = ~J~T7ncr 17 ’ 2- 1 

c (e -1) 

where : 

h is Planck’s constant = 6.625 x 10 34 Ws 2 , 
c is the speed of light, 

k is Boltzmann’s constant = 1.380 x 10 23 J-K" 1 ’ 
v is radiation frequency, Hz, and 
T is temperature in Kelvin. 

Mbb ( v ,t) = Spectral Radiant Exitance, also called Radiant Emittance, is defined as 
the power per unit area per unit wavelength interval at a particular wavelength [Ref.8]. 
This equation can be simplified using defined constants and expressed in terms of 
wavelength, X: 
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2.2 



M 



C, 



BBW A\e c ' UT 



1 ) 



Where : 

Ci = 2nhc 2 = 3.741 x 10 4 W cm' 2 Jim 4 ' 

C 2 = ch/k = 1.438 x 10 4 Jim K, and 
X = is the wavelength, Jim. 

Both equations express the spectral radiant exitance/emittance of a source 
radiating in a vacuum. We see in Figure 2.2 the spectral radiant emittance of a blackbody 
at different wavelengths. 




Figure 2.2 Radiant Emittance of a Blackbody [Ref.8] 
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When we are dealing with non-blackbodies, we have to consider an additional 
factor, which introduces the special characteristics of the material composition and 
surface condition. It is called “the emissivity constant” and represents the ratio between 
the emitted radiation of any given temperature radiator and the emitted radiation of a 
blackbody of the same temperature. 



Mb = £ Mbb , 



2.3 



where : 

£ = is the emissivity of the non-blackbody. 

2. Stefan-Boltzmann Law 

The total radiated energy from a unit area of a surface per unit time over a 2 K 
steradian solid angle is called the total radiant exitance of that surface. The dependence of 
total radiant exitance on the absolute temperature of a blackbody is described by Stefan- 
Boltzmann’ s law: 

M ( t) = a T 4 2.4 

where a is the Stefan-Boltzmann constant = 5.67 x 10' 12 W-cm' 2 K' 4 

The same relationship as before has to be considered if the object is a non- 
blackbody. So, radiation falling on a body is absorbed, reflected, or transmitted and 
because of the law of conservation of energy, 
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CX+ P + T = 1 



where a is the absorptance, p is the reflectance, and x is the transmittance with all of 
them wavelength-dependent. 

3. KirchhofPs Law 



Kirchhoff observed that, at given temperature, the ratio of a body’s radiant 
exitance to its absorptance is constant, and equal to the radiant exitance of a blackbody, at 
the same temperature: 

— = W 0 This is Kirchhoff s law. 

a 

Therefore, if a body is a good absorber of radiation, its emissive capability, as 
expressed by W, must be good too, since the ratio is constant. That is why a good 
absorber is also a good emitter. 

Thus, 



£QT A 

a 



= oT 4 and, therefore, 



2.5 



e = a, the emissivity of a body is equal to its absorptance and 
depends on: 



* the nature of the body. 

* its temperature. 

* its surface finish. 

* the wavelength. 
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Another consideration is the fact that we often work with “atmospheric windows,” 
which means that we can subdivide our spectrum into specific sections of interest. In 
these cases we can estimate the “In-band Fractional Power” or in-band radiance “N” 
[Ref.5], using the “Universal Blackbody Curve” or the proportional radiation table, with 



N = 



£(JT 4 



71 



-A<?, 



2.6 



Where Aq = difference of two q values for the limiting wavelengths from the 
proportional radiation table (See Appendix A). The proportional radiation q is a function 
of the wavelength and temperature. It represents the fraction of the radiant exitance 
emitted by a blackbody at temperature T at all wavelengths up to the specific value of 
wavelength, 

4. Lambert-Beer Law 

In order to “see” the target, we have to know or detect the radiance seen by the 
sensor, after propagation through the atmosphere. In order to do this, we have to know 
how the atmosphere interacts with the optical ray on its way to the sensor. It has been 
stated that the transmittance is a strong function of the wavelength, which can be 
described by the Lambert-Beer law: 

T( x ) — e , 2.7 
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where, 



T( x ) = Spectral atmospheric transmittance, 

p(t) = extinction coefficient, 

R = Path length, 

The extinction coefficient depends on the atmospheric molecular composition and 
aerosol concentration. It is very sensitive to the specific frequencies of the molecular 
vibration and rotation transitions. Of special interest for IR emission are energy 
transitions of the vibrational type that produce spectra in the band from 0.3 to 30 pm 
[Ref. 1 : p.190]. 

The infrared radiation will be attenuated as it propagates through the atmosphere 
by the processes of absorption and scattering. The attenuation is characterized by the total 
extinction coefficient, 

P = Pa + Ps. 2.8 

where, 

p = total extinction coefficient, 

p a = extinction coefficient for total absorption, 

p s = extinction coefficient for total scattering. 

Each of these two extinction coefficients has two other components; one from the 
molecules of the air and the other from aerosol particles suspended in it. Scattering by 
aerosols and molecules has a large effect in the visible region, while absorption 
dominates in the infrared region of the electromagnetic spectrum. As the EO systems are 
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designed to operate mainly between 8 and 12 pm, it is important to predict the 



transmittance of the atmosphere as a function of wavelength and weather conditions. 
These two extinction coefficients are described by the equations: 



|^a — K m + K a y 


2.9 


|^s = Gm + G a , 


2.10 



where, 

K m = molecular absorption coefficient, 

Kz = aerosol absorption coefficient, 

G m = molecular scattering coefficient, 

G a = aerosol scattering coefficient. 

Earth’s atmosphere is a mixture of many gases with varying characteristics of 
absorption, emission, and scattering of optical radiation. Each of these gaseous quantities 
varies with altitude, time, and space as a function of geographical region. [Ref. 4: p.71] 
So, the relative values of the four coefficients will change depending upon the region 
where the measurement is made. 

Due to the large number of parameters involved in optical transmission a range of 
computer models have been developed for computation of atmospheric transmittance and 
atmospheric path radiance for various defined meteorological scenarios. These models 
include LOWTRAN, MODTRAN, HITRAN and variations of those. 
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Gathering all these concepts and theoretical and experimental laws we can get an 
idea how detection scenario looks in the 1R arena. Broadly speaking we have three main 
components: a target and its surroundings, a sensor system on a moving platform, and an 
atmospheric path between them. 

According to Holst, the 1975 NVL and FLIR92 models are the main analytic 
tools for deriving system requirements and predicting performance for FL1R systems. The 
1975 NVL model was developed for serial and parallel scanning thermal imaging systems 
that existed in the 1970s. FLER92 is built upon the framework of the 1975 NVL model. 
The three-dimensional noise model and FLIR92 are in symmetrical formats. This permits 
a format that is easier to understand than the 1975 NVL model, although both provide 
similar results [Ref. 2:p. 389-391]. 



C. FLIR SYSTEM PERFORMANCE 

The purpose of a FLIR system is to detect the arrival of IR energy, electronically 
process the detected energy, and display the results to the FLIR operator. There are two 
standard parameters for describing FLIR performance: Minimum Detectable 

Temperature Difference (MDTD or MDT) and Minimum Resolvable Temperature 
Difference (MRTD or MRT). Each of these provides a sensitivity measurement of the 
FLIR system. 
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1 . 



MRTD 



The MRTD method of FLIR system description is the most commonly used. The 
definition of MRTD is “ the smallest temperature difference in a standard periodic test 
pattern for a given spatial frequency that is resolvable to an observer over an unlimited 
viewing time" [Ref.6:p.242]. So, it gives the AT between the target and ambient 
temperature of a standard 4 bar (7:1 aspect ratio) chart required to make the bars just 
resolvable as a function of the spatial frequency. 



Shumaker gives a typical form for the MRTD equation. [Ref. 7: p.8-52] 



MRT(v) = 



2 * SNRT * NET * 

MTF(v) 



Yx 



Y * A * A, 



Yi 

[t e *F r *N 0 S *NjX, 



2.11 



where, 

SNRT = signal to noise ratio threshold 
NET = Noise Equivalent Temperature 

MTF = Modulation Transfer Function 
N = Overscan ratio 

OS 

N ss = Serial scan ratio 

L = Length-to- width ratio of the bar 

v = Spatial frequency 

t e = eye integration time in seconds 

F r = Frame rate 
p x = Noise filter factor 
A x = in-scan detector subtense 
A v = cross-scan detector subtense 
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For calculations purposes it is convenient to express MRTD as a function 



of spatial frequency f in a parametric form. 



MRTD = SL + SC*2* 




exp 



(* 



nf.R'f 



2.12 



where SL is the sensitivity limit, SC is the sensitivity contrast and ER is the Equivalent 
Resolution. All three constants are derived for any given FLIR by fitting the equation 2.12 
to the measured or calculated MRTD of the system for a 4 bar target pattern with aspect 
ratio equal to 7. FLIR constants are not generally available for specific military FLIR 
instruments but can be deduced from suitable generic models. Since MRTD is a function 
of spatial frequency, the maximum range for classification or recognition is the range at 
which the ship/background AT matches MRTD at the spatial frequency given by the 
appropriate criterion (e.g., the Johnson Criterion). 

2. MDTD 

The MDTD definition of a FLIR system is similar to the MRTD, but the target is 
represented as a large square and MDTD is not as dependent on the eye parameters. The 
definition of MDTD is the minimum detectable temperature difference between a large 
square target and its background that just makes the square target detectable. Shumaker 
also provides a typical form for MDTD [Ref. 7: p.8-55]. 

The calculations of MDTD for a generic FLIR system can be described by a curve 
fit equation as shown in equation 2.13, 
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2.13 



MDTD = SL + 



SC 
TAS 2 



(fAS 2 +ER 2 




where SL, SC and ER are empirical FLIR constants as before. TAS represents the 



target angular subtense in mRad, which equals 10 6 




Aj is the projected area of the 



target, and R is the range. 

The effective temperature difference, AT, is the temperature difference required of 
two blackbodies to produce the actual in-band average radiance difference between a 
target and its background. Shumaker provides a detailed discussion of all the factors that 
must be considered in determining the nature and magnitude of thermal signatures, and 
relates them to a AT that can be used to characterize targets for analysis [Ref.7:p.2-19]. 



3. Johnson Criteria 

In 1958, Johnson conducted an experiment using eight military vehicles and one 
observer to develop a relationship between spatial frequency and a target. During the 
experiment, Johnson placed tri-bar patterns next to military targets as he varied 
conditions. He noted what tri-bar frequency could just be resolved when a given visual 
discrimination task was accomplished. Then, he related the visual discrimination levels to 
the Air Force tri-bar target frequencies using the minimum dimension of the target as a 
reference. From this experiment, and some improvements to relate the results to FLIR 
systems the Johnson Criteria was established. The table 2.1 contains the two dimensional 
cycle criteria for a 50% probability of task discrimination that are in current use. 
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Discrimination Level 


Definition 


2-D Cycles 


DETECTION 


An object is present (Object vs. Noise) 


0.75 


CLASSIFICATION 


Distinguish targets from similar sized non-targets 


1.5 


RECOGNITION 


The class to which the object belongs (tank.Truck, 
Destroyer) 


3.0 


IDENTIFICATION 


Determines the specific military designation (T72 
tank.DD963, CVN69) 


6.0 



Table 2.1 Two-dimensional Johnson Cycle Criteria. 



In order to establish the bar-chart equivalent of a specific discrimination task, 
selecting the proper target size is critical in obtaining meaningful results. In Johnson’s 
original work, he chose the minimum target dimension as the target’s “ characteristic 
dimension Lawson and Moser found that using the average target dimension (square 
root of its area) as the “ characteristic dimension" led to consistent predictions for targets 
regardless of aspect, while use of the minimum dimension led to erroneous conclusions. 
The area of the target that is of concern to the FLDR analyst is that projected area normal 
to the line-of-sight from the sensor to the target. This area is dependent on the target 
dimensions, the observing aspect, and any obscuration of the target. 

Shumaker provides detailed descriptions and physical dimensions for a variety of 
military targets [Ref.7:p.2-8/2-14]. 
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III. COMPETING MODELS 



A. TACTICAL DECISION AIDS 

The Tactical Decision Aids (TDA) are special tools designed for assisting a 
decision maker by assimilation and convenient presentation of data and analysis of 
tactical problems beyond what is feasible by humans in timely fashion. Since the 
advances of the technology go on and on, the weapon systems have reached a high level 
of development that sometimes makes them very difficult to use in the operational 
environment. This is one of the main reasons why we have TDAs, [Ref.9: p.1-1]. 

According to Holst, Electro-optical imaging system analysis is a mathematical 
construct that provides an optimum design through appropriate tradeoff analyses. A 
comprehensive model must include the target, the background, the properties of the 
intervening atmosphere, the optical system, detector, electronics, and the human 
interpretation of the displayed information. While any of the components can be studied 
in detail separately, the whole Electro-optical imaging system cannot. [Ref. 2 : p.12-13] 

An EOTDA typically consists of input and output processors and three main 
models: a Target Model that determines the inherent signal that emanates from scene 
objects, an Atmospheric Model which estimates the degradation of signal from target to 
sensor due to weather conditions, and a Sensor Model which determines the range at 
which an object can be detected or at which a weapon system can lock-on. This thesis 
describes each one of these two codes (under study) based upon these three main or 
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fundamental parts plus the way that each one presents the outputs. See Figure 3.1 on page 
58 for a Functional Diagram of a general TDA. 

B. ACQUIRE MODEL 

ACQUIRE is an analytical model that predicts target detection and discrimination 
range performance for systems that image in the visible, near infrared and infrared 
spectral bands. Broadly speaking it is a FORTRAN program designed to run on IBM- 
compatible PCs and Unix workstations. Also, it is intended for experienced system 
analysts who are knowledgeable of IR technology and the basic methodology applied in 
the model. 

This thesis uses ACQUIRE version 1 dated May 1995. The model was developed 
by the U.S. Army Night Vision and Electronic Sensors Directorate, at Fort Belvoir, 
Virginia. As for other Electro-optic system range performance models, three main 
components are built in: 

1. The Target Model 

The program calculates the probability-range relationship based on input target 
characteristics (length, height, and target-background temperature difference), 
atmospheric conditions, MRTD curves, and a 2D formulation of the Johnson’s Cycle 
Criteria methodology. Basically, the target is represented by a temperature difference 
against background. The target signature is either the target contrast (visible or near IR 
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imager) or “delta T” the zero-range temperature difference between the target and its 
background (thermal imager). 

The model includes range prediction methodologies for both target discrimination 
and target spot detection. For a target imbedded in a nonuniform or cluttered background, 
discrimination of characteristics that separate the target from the background is required 
in order to perform detection. For a target against a uniform background, detection may 
occur when the signal-to-noise ratio (SNR) on the display element representing the target 
exceeds the SNR of the background. 

Target discrimination is usually divided into three levels: detection, recognition, 
and identification. Detection can be defined as the determination that a target of potential 
military interest is present within the sensor field-of-view. Recognition would be the 
discrimination between specific objects within a class of similar objects (truck, tank, 
APC, etc). Identification is defined as the discrimination between specific targets (T-72, 
Ml, etc) [Ref. 10:p.l], 

Due to advances in technology and improvements to FLIR systems, Johnson’s 
Criteria have had to be updated. A series of perception tests conducted by the U.S. Army 
CECOM Center for Night Vision and Electro-optics Visionics Division demonstrated that 
both horizontal and vertical resolution are equally important to determine the ability of 
test subjects to perform visual discrimination tasks. This correction was made to the 
Johnson Criteria, which generated a new level of target discrimination, classification that 
can be defined as discrimination between general classes of vehicles [Ref. 10: p.13]. 



21 



ACQUIRE User’s Guide defines the probability of target discrimination as a 



function of the number of equivalent cycles resolved on the target by the sensor. For a 
given target at specific range and apparent signature, the number of cycles resolved is 
determined by either Minimum Resolvable Temperature Difference (MRTD or MRT) or 
Minimum Resolvable Contrast (MRC), and is given by: 



C D is the Characteristic Size of the target (Critical Dimension) in meters, R is the range 
to the target in kilometers and f R is the frequency in cycles per miliradian resolved by the 
sensor. ACQUIRE uses an empirical curve fit to the Target Transfer Probability Function 
(TTPF) that is given by: 



where E = 2.7 + 0.7 ( n/n 50 ). 

\n/n 50 ) is the number of cycles required to be resolved in order to achieve a 50% 

probability of target discrimination (detection, classification, recognition, or 
identification) [Ref. 9]. 

In the case when a target is against a uniform background, the detection range is 
calculated using either Minimum Detectable Temperature Difference (MDTD) or 
Minimum Detectable Contrast (MDC) and threshold signal-to-noise ratio. 




3.1 



where, 




(n/n 50 ) E 



3.2 
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ACQUIRE includes 23 standard targets (small land targets) to select from the 
target_lookup table, if the target is not included, a target data group must be created. This 
data group is used to specify the size, as its critical dimension, and signature of the target. 

2. The Atmospheric Model 

To reach the IR detector, the radiant flux from the target has to pass through the 
atmosphere. The earth's atmosphere is filled with absorbing agents, components of the 
atmosphere, which extract energy from the supply of radiation (attenuation). Also, the 
small particles suspended in the atmosphere cause scattering and emission of IR radiation 
in all directions. 

The most significant absorbers in the atmosphere are water, carbon dioxide, and 
ozone. On the other hand, scattering is produced by molecules of the air and aerosol 
particles suspended on it. The ratio of the intensity of usable radiation passing through a 
body of air and entering an optical system to the intensity originating from the target is 
known as the atmospheric transmittance [Ref. 7 :p.3-l]. 

ACQUIRE has two options for modeling atmospheric transmittance. The first 

option is using Beer’s Law, in which the atmospheric transmittance (T a ) is given by: 

X 3.3 

where |l is the attenuation (extinction) coefficient (km -1 ) and R is the range (km). The 
extinction coefficient includes the molecular scattering coefficient, the molecular 
absorption coefficient, the aerosol absorption coefficient, and the aerosol scattering 
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coefficient. For a finite bandwidth, Beer’s Law with a single constant extinction 
coefficient frequently does not adequately represent the extintion. 

The second option, recommended in the ACQUIRE Users Guide, is to specify 
atmospheric transmittance as a function of range. Several models exist for computing the 
transmittance; the most commonly used is LOWTRAN. Currently, Lowtran is being 
retired in favor of MODTRAN. This method requires the standard meteorology data set. 

The effects of smokes and obscurants are included by specifying a mass extinction 
coefficient and a concentration length. ACQUIRE, then modifies (attenuates) the target 
signature due the atmospheric transmittance and smoke transmittance. 

3. The Sensor Model 

The sensor model in this code is selected according to its specific performance 
parameters. They are selected as MRT, MDT, MRC, or MDC from a look-up table. Each 
sensor has a specific character string ED which is required as input data. The sensor IDs 
must match with sensor performance curve data that is stored in a special data file. 
ACQUIRE also lets users define specific sensors when they are not included in the 
lookup table. In this case, the horizontal field of view is required in order to calculate 
search performance; it must match the input MRT or MRC data points. 

C. WIN-EOTDA 

The Naval Research Laboratory, in Monterey, CA, developed Win-EOTDA 
version 1.0. This product was derived from the Mark HI Electro-Optical Tactical Decision 
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Aid (EOTDA), which was originally developed by the USAF Phillips Laboratory. 
This thesis uses Win-EOTDA version 0.7 dated June 1997, which is a special version just 
for experimental purposes [WinEOTDA Program Help]. 

Windows Electro-Optical Tactical Decision Aid (WinEOTDA) is a computer 
model that predicts the performance of Electro-optical systems, based on environmental 
and tactical information. The model consists of three microcomputer-based programs 
supporting infrared (8-12 |im), visible (0.4-0.9 |im), and laser (1.06 p.m) systems. Each 
program is comprised of three sub-models: an atmospheric transmission model, a target 
contrast model, and a sensor performance model. 

WinEOTDA has been designed to simplify the process of predicting performance 
of Electro-optical weapon systems and night vision goggles (NVG). The input and the 
outputs are presented on a single easy to read form, shown as Figure 3.2. It has a main 




Figure 3.2 WinEOTDA Main Screen 
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screen that contains both input and output in a graphical representation, both of them 
available at a single click of the mouse. 

Input data can be changed from the main menu, the toolbar, and by clicking on hot 
spots on the main screen. Data entry forms will appear and guide users through the input 
process. Output can be found on the main screen and on a separate form that appears 
when you select “Output Form” from the main menu or toolbar. Users will navigate 
among the various screens in the model by using both drop down menus and navigation 
buttons, just like any other Windows application. 

The WinEOTDA main screen leads to multiple forms (and sub-models) that allow 
users to access a user-friendly environmental tool for performing range predictions in a 
variety of scenarios. 

1. The Target Model 

In the target model, the program calculates the difference in radiance between 
target and background over the wavelength band from 8 to 12 p.m. This radiance 
difference is then converted to a temperature difference. The thermal model converts the 
physical temperature of the target and background to an equivalent blackbody 
temperature. 

The thermal model is based on the Target Contrast Model #2 (TCM2), which 
treats the target as a three-dimensional network of nodes that exchange heat with one 
another and with the environment, providing a detailed thermal signature. The number of 
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nodes per target varies from 17 nodes for the POL tank to 68 nodes for the Apache 
Helicopter (AH-64). 

The principal phenomena interacting to produce the thermal scene are radiative 
heating and cooling; mass and heat transfer effects of evaporation, condensation, 
sublimation, and precipitation; the thermal properties of the target and background; and 
the thermal transfer by wind [Ref. 1 l:p. A-l], 

The TCM2 computes a mean target temperature and it also identifies the hottest 
and coldest of the visible nodes. It then identifies which of these has the greater contrast 
with the background. This identified node gives the target contrast temperature 
difference. 

The WinEOTDA includes 18 standard military (tanks, frigate, gunboat, helicopter, 
power plants, etc) targets to select from, as well as eight classes of generic targets, which 
the user can specify. Each target is associated with up to eight types of backgrounds. The 
backgrounds available are vegetation, soil, snow, water, concrete, asphalt, swamp, and 
rocky field. Each one of these backgrounds is further described by moisture, coverage, 
depth (for water background), and density of the surface type. The WinEOTDA is usually 
used with multiple background instances (up to three) for each run. The program serves to 
use the first entered background as the primary background, which is used to compute 
reflected ground radiation. Each entered background is not considered independently. It is 
important, therefore, to enter the most representative background first. 
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2 . 



The Atmospheric Model 



The atmospheric model evaluates the attenuation of the signal by atmospheric 
constituents between target and sensor. The attenuation consists of four components: 
molecular, aerosol, precipitation, and battlefield-induced contaminants. 

Molecular attenuation is due primarily to absorption by water vapor and is related 
to air temperature and moisture. It effects directly the Mid IR and Far IR bands as is 
shown in the following Table: [Ref. 13 :p.7-3]. 



WAVELENGTH REGION 


ATTENUATION COEFICIENTS 


Visible (0.4 - 0.7 pm) 


Aerosol scattering - Molecular scattering 


Near IR (0.7 - 2.0 pm) 


Aerosol scattering - (Molecular absorption only for 
specific wavelength) 


MidIR (3.0 -5.0 pm) 


Molecular absorption (water vapor) - Aerosol scattering 


Far IR (8.0 -12 pm) 


Molecular absorption(water vapor) - Aerosol scattering 



Table 3. 1 Importance of attenuation coefficients according to wavelength. 



The combined transmittance resulting from the four types of attenuators is 
accurately evaluated by the model for a range of four kilometers. Transmittance at other 
ranges is estimated by assuming the exponential extinction law (Beer’s law). 

The WinEOTDA model contains 17 different types of aerosols from two basic 
models- LOWTRAN versions 5,6, and 7 use nine aerosol models and eight more from the 
Navy Aerosol size distribution model (NAM). The following Table shows the nine 
aerosols models contained in WinEOTDA [Ref. 10: p. 127]. 

Most of the aerosol model is accessed by using the entering aerosol model data. 
The rural model describes the boundary-layer background aerosol found in continental air 
masses (microscopic soil particles and manmade processes). The urban model describes 
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RURAL 



URBAN 

MARITIME 

TROPOSPHERIC 

DESERT 

NAVY MARITIME 
ADVECTIVE FOG 
RADIATIVE FOG 
CAMOUFLAGE SMOKES 
Table 3.2 Aerosol Models 

the concentration of particles produced in urban and industrial complexes. The 
tropospheric model characterizes aerosols found in very clean air masses and in the free 
atmosphere. The desert aerosol model characterizes aerosols found in the boundary layer 
of desert, arid, or semi arid climatic regions. 

The Navy Maritime model describes the aerosol found in oceanic environments. 
The main difference from the other maritime model is its wind speed dependence. There 
is a sub menu with nine options that depend on air mass history and the average wind 
speed over the target area during the 24 hours previous to the aerosol model forecast time. 

3. The Sensor Model 

The sensor model describes the performance of the sensor in terms of the 
minimum resolvable or minimum detectable temperature difference as a function of 
spatial frequency. A user defined FLIR model can be constructed using the physical 
parameters of the system: overall field of view FOV, instantaneous field of view IFOV 
along and across scan, detector sensitivity, NET, system modulation transfer function 
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MTF, etc, from which MRT and MDT are deduced as function of spatial frequency 
(cycles/mrad). 

The range performance model then determines the range at which the target 
contrast temperature differences modified by atmospheric transmission equals the sensor 
MRT or MDT. The relation between target spectral frequency (for given task of detection, 
recognition, etc) is defined by target size and the Johnson criterion. 

The sensor performance model determines the range at which the actual contrast 
(temperature) signal received by the sensor equals the threshold signal for detection or 
lock-on. Normally, the signal decreases as the target-to-sensor range increases. The 
threshold signal is determined by the target’s apparent size (angular subtense) as viewed 
from the sensor. The angular subtense is equal to the critical dimension of the target 
divided by range to sensor. The smaller the angular subtense, the larger the threshold 
value of the signal. 

WinEOTDA supports three sensor types (IR, TV, or Laser). This thesis only treats 
IR sensors. A number of sensor data files are supplied with the program. For some of the 
supported sensors, system parameters and/or predicted ranges are classified when 
associated with the sensor name. Each sensor has a unique three-digit number associated 
with it. Sensor IDs are divided into two categories: standard IDs reserved for sensors 
supplied with the program and additional IDs for user-defined sensors. The list of 
software-supplied sensors and IDs is in a confidential appendix of the User’s Manual. 
Neither the list nor the confidential appendix was available for the author of this thesis 
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because of the classification level. This thesis utilized the generic sensor model genl and 



gen2 included in WinEOTDA and one user-defined model NACIT9. 



D. SCENARIO 

Since EOTDA assess the performance of air-to-ground weapon systems and 
ACQUIRE does the same task for small land targets, particular scenario conditions were 
selected for achieving and setting the comparison between these two codes. The scenario 
will be described according to the main sections reviewed previously. 

A land target represented by a Soviet T-72 tank and a Gunboat as an air-to-sea 
target were selected based on the following criteria: 

• Target related with TDA Code designed for. 

• Target in a lookup table for both TDA codes. 

• Target data availability in the references. 

Operational sensor performance data were not available for use in WinEOTDA 
and ACQUIRE due to classification issues. Performance was estimated based on 
published physical descriptions of the AN/AAS-36 IRDS system. A general description 
of this sensor can be found in the thesis by Kreitz [Ref. 1 l:p. 17-21], The NACIT9 MRT 
model was generated using these data. The two generic models genl and gen2 were also 
used. 



31 



The atmospheric data were selected according to the typical operational scenario 
for each target selected. Therefore, a desert atmospheric data set for a T-72 tank target 
obtained from Strike Warfare Weapons School at Fallon Naval Air Station, Nevada 
[Ref. 16] was used. Special open maritime atmospheric conditions for a Gunboat target 
were selected from the database of the PREOS92 experiment at Monterey Bay, performed 
by Naval Command, Control and Ocean Surveillance Center (NCCOSC), Research, 
Development, Test and Evaluation Division (NRaD) and Naval Postgraduate School 
(NPS) [Ref. 13], 

The following Table summarizes the data used as input for these targets in the 
ACQUIRE and WinEOTDA codes. The selection of ACQUIRE with target T-72 and 
WinEOTDA with target Gunboat is just an example. Each target was exposed to each 



code with similar parameters and under similar conditions. 



DESCRIPTION 


ACQUIRE: T-72 


WinEOTDA: GUNBOAT 


TYPE OF TARGET 


LAND TARGET 


AIR-TO-SEA TARGET 


LOCATION 


39° 14’(N) - 1 18° 14’ (W) 


36° 47’ (N)- 122° 18’ (W) 


DIMENSION (mts) 


H:2.16 L:9.33 W:3.35 F:64 


H:8.8 L:41.5 W:9.75 F:60 


AEROSOL 


DESERT MODEL 


NAVY MARITIME MODEL 


SENSOR SYSTEM 


AN/AAS-36 (IR) 


AN/AAS-36 (IR) 


TRANSMITTANCE 


Calculated by Lowtran 6 


Not required 


TEMP(Min/Max/Min) 


29/47/29 F 


42/55/46 F 


Table 3.3 Summary of data input for both coc 


es: ACQUIRE and WinEOTDA 



AN/AAS-36 FLIR represented by NACIT9 model 
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IV. DATA COMPUTATIONS AND RESULTS 



A. GENERAL 

ACQUIRE and WinEOTDA have their particular ways to handle data inputs, 
computing or calculating algorithms and the output data. This work is not intended to 
describe in detail the way that each code does this; it would be better to have a Technical 
Manual on hand for this purpose. Instead of this, the main modules of each code are 
presented with the data processing used to get the predicted target detection range for 
each for equivalent operator inputs. It is necessary to point out that both codes can 
perform range prediction in three well-defined bands; however this work just deals with 
the infrared band (8 to 12 pm). 

According to Holst, performance of an optical system is optimized by performing 
a series of tradeoff analyses as a function of a variety of variables. As we saw previously, 
typical variables include field of view, atmospheric transmittance, target size, target 
intensity, and stabilization. Each tradeoff provides a different view of overall 
optimization [Ref. 2:p 441], These tradeoffs will be pointed out in this section. 

B. ACQUIRE INPUTS 

The Army code ACQUIRE requires a data input file with the target, sensor and 
scenario parameters, which must be created before running the code. This data file must 
be an ASCII text file with a specified format (See Appendix B). As was pointed out in 
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chapter two, the target can be defined using a targetjookup table or a user defined data 
group. The two methods were used in this work. 

The target_lookup entry format used was: 

>target_lookup 

target_id T72 

aspect F F_or_S 

signature 1.25 degrees_C_or_contrast 

The target_id selects a T72 target from the lookup table. There are 23 small land 
targets in this table; also, there is a default value for a T62 target. The aspect represents 
the position of the selected target; it is related to the critical dimension directly. For the 
Front (F) aspect the critical dimension value is 2.55 meters and for the S aspect the 
critical dimension is 3.60 meters; these values are set in the database for each target. The 
code was run twice, once in each aspect to discriminate the difference in the performance. 
The results are included in the Appendix B. The signature is the difference in 
temperature between the target and the surrounding background. It was set as the default 
value for a thermal imaging system (1.25 degrees centigrade). 

The other selected target was a Gunboat. In this case it was necessary to define the 
target characteristics because it is not represented in the internal lookup table. The format 
used was: 

>target 

characteristic_size 10.12 meters 

signature 4.6 degrees _C 

The Gunboat is then represented by these two parameters. The 

characteristic_size represents the critical dimension of the target as defined previously. It 
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was computed for both aspects Front and Side using the following equation from 
Shumaker [Ref.7 :p.2- 15]: 

A t = Iw sin 6 cos <j> + hw cos 6 cos 0 + hi cos 6 sin (j) 4.1 

where, 

~\ 

/ is the target length 

iv is the target width 

h is the unobscured target height > of an equivalent Parallelepiped 
$ is the azimuth angle 

6 is the elevation angle 

The following Table shows the parameter values of the Gunboat target based on 
the dimensions of the Research Vessel POINT SUR (R/V POINT SUR) which was used 
as a platform for the PREOS experiment in 1992. The Table also includes the critical 
dimension values computed by Equation 4.1. 



Ship length 1 


41.5 meters 


Ship width w 


9.75 meters 


Ship height h 


8.80 meters 


Elevation angle 6 


2.80° sin -1 (altitude / range) 


Azimuth Angle Side Aspect (f> 


0° 


Azimuth Angle Front Aspect <j> 


90° 


Projected Area A T Side Aspect 


384.16 m 2 


Projected Area A T Front Aspect 


102.41 m 2 


Critical Dimension Side Aspect 


19.60 meters 


Critical Dimension Front Aspect 


10.12 meters 



Table 4.1 Critical Dimension of Gunboat target. 
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The signature parameter is the temperature difference between the target and its 
surrounding background; it was set to 4.6 degrees as the average temperature of the ship 
and the background during the selected day on August 4,1992. [Ref. 12:p.57-59] 

The environmental effects are considered in ACQUIRE by modeling the 
atmospheric transmittance and the effects of smokes and obscurants. The ACQUIRE 
User’s Guide suggests getting the atmospheric transmittance as a function of Range. 
Using LOWTRAN 6, an atmospheric propagation model, a set of values was obtained 
relating to the specific environment where those targets were supposed to operate; a 
desert environment and a maritime environment. Table 4.7 and 4.8 (pag.59-60) shows the 
parameters used as input in LOWTRAN 6. This is the only way in which ACQUIRE 
introduces meteorological and regional data into calculations. The format used was: 
>band-averaged_atmosphere 

#_points: 20 Km transmittance 

0.000e+00 1.000e+00 

The transmittance values are included as part of the written program for each 
target, and they are shown in Appendix B. Smokes and obscurants were not considered. 

Another main section in the ACQUIRE input program is the sensor. It is 

considered as the vital part of the program. For this thesis several programs were written 

and run using the generic sensors (genl, gen2) included with the ACQUIRE model. Also, 

a sensor data file was developed by the author from the available data of the AN/AAS-36 

and theoretical values. The required input data defining a sensor are the MRTD and 

MDTD values as functions of spatial frequency, expressed in terms of empirical constants 
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equivalent to SL, SC, and ER of equation 2.12. An appropriate set of constants was 
obtained by inserting the available physical parameters of the AN/AAS-36 FLIR in 
equation 2.11 with theoretical estimates for the terms not directly available. The 
corresponding set of FLIR parameters SL, SC, and ER was determined by matching 
equation 2.12 to equation 2.11 empirically for the useful range of spatial frequency. The 
optimal values of the parameters were determined iteratively using an Excel spread sheet 
to match the outputs of equation 2.12 to equation 2.1 1. Table 4.2 shows the set of values 
of MRTD as a function of spatial frequency, which was developed to be used as sensor 
inputs for each code. The estimated MRTD curve and MDTD curve for calculating range 
predictions were derived from Shumaker’s example 8-13 and 8-16 [Ref. 7 ] for a generic 
FLIR system in WFOV. The detector angular subtense Ax and Ay were chosen as 0.45 



mRad to represent a modeled FLIR system in WFOV. The Noise Equivalent Temperature 
Difference (NET) was chosen as 0.140 K, between 0.1K (Shumaker example) and 
0.250K (Moser’s reference). Figure 4.1 is the Excel data sheet including the graph of 
MRT curves on page 61. 





From Shumaker 


Result 




0.25mFtad 


ImRad 


0.45mRad 


MRT= 


0.010018 


0.057341 


0.017026 




0.012722 


0.120642 


0.024775 




0.015714 


0.268442 


0.034768 




0.019114 


0.675595 


0.048469 




0.023063 


1.978021 


0.068144 




0.027737 


6.835714 


0.097468 




0.033359 


28.12041 


0.142616 




0.040214 


138.4451 


0.214273 




0.054719 


2599.65 


0.436981 




0.133215 


56697152 


4.230673 




0.255343 


2.07E+1 1 


25.26968 



Table 4.2 MRTD as function of spatial frequency for “generic” FLIR system 
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The sensor data file contains the data (MRT and MDT) as a function of the spatial 

frequency for gen 1 and gen2 from the ACQUIRE database (included with the code) to be 

used for computing the prediction ranges. The sensor file is enclosed as a part of 

Appendix B. The format used was: 

>sensor_lookup 

data_file_name sensor.dat 

sensor_id genl 

performance_mode MRT MRT_MDT_MRC_or_MDC 

Using this format, the code was run three times, first one with the data from genl 

and second with the data from gen2. The results show slight differences. The third run 
was with the prototype sensor data (calculated from AN/AA-36 and theoretical 
parameters) which was included in the same sensor data file under the name “NACIT9”. 
These runs used the same format as above, but with the sensor_id changed from genl to 
gen2 and NACIT9. 

The program files with all needed parameters were arbitrarily designated as 
NACm, NACIT11, NACIT2, NACmi, NACIT3, and NACIT31. All of them are 
included in the Appendix B. ACQUIRE generates an output file with the same name as 
the input file plus an extension. 

C. ACQUIRE OUTPUTS 

From the input files reviewed previously, ACQUIRE generates an output file for 
each one. Each output file has an extension - rl or r2 - The rl file is used to verify the 
performance with respect to a specific probability requirement. The r2.file lists in tabular 
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form the discrimination probability given target range, which is the way that ACQUIRE 
calculates performance. Six output files were obtained after the ACQUIRE computations, 
and these are included in the Appendix B. 

The NACITl.rl and the similar files are basically divided into three sections. The 
first section is a summary of the input and default values used by ACQUIRE. The second 
one is a detailed examination of the intermediate result files and is presented as text 
messages. Figure 4.2 shows an example of this section, for the Gunboat target at Side 



aspect. 

MESSAGES 

2D MRTD from lookup table 



cyles/mrad 


MRTD 


2 . 970e-001 


1 . 02 0e-002 


4 . 3 00e-001 


1 . 460e-002 


5 . 990e-001 


2 . 100e-002 


8 . 200e-001 


3 . 02 0e-002 


1 . 091e+000 


4 . 340e-002 


1 . 3 90e + 000 


6 . 23 0e-002 


1 . 692e+000 


8 . 960e-002 


1. 987e+000 


1 . 28 8e-001 


2 . 265e+000 


1 . 851e-001 


2 . 522e+000 


2 . 660e-001 


2 . 7 54e+000 


3 . 823e-001 


2 . 959e+000 


5 . 494e-001 


3 . 148e+000 


7 . 896e-001 


3 . 300e+000 


1 . 13 5e+000 


3 . 440e+000 


1 . 63 le+000 


3 . 567e+000 


2 . 344e+000 


3 . 649e+000 


3 . 369e+000 


3 . 723e+000 


4 . 841e+000 


3 . 77 8e+000 


6 . 958e + 000 


3 . 834e+000 


1 . 000e + 001 



Sky-to-ground ratio defaulted to 1 for thermal systems. 

Last range for input atmospheric transmittance data is less 
than maximum range . 

Spline fit returned a negative atmospheric transmittance; add 
more points to input transmittance data array. 

Loadline does not intersect MRTD/MRC above range 15 km, extend 
curve to lower frequencies and temperatures. 

SENSOR 

genl from sensor.dat 

TARGET 

characteristic dimension: 19.6 meters 
inherent signature: 4.6 degrees C 



Figure 4.2 ACQUIRE typical output (Messages) 
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The third section is the tabulated result of the run, called “Observer Ensemble 



Performance Range Given Probability...”. The probabilities are given in five- percent 



steps (5%) from 0.95 percent down to 0.05 percent (5%). The following Figure shows a 



typical output file that is generated from ACQUIRE, 



OBSERVER ENSEMBLE PERFORMANCE 



RANGE GIVEN 


PROBABILITY. . 












WFOV 


NFOV 


NFOV 


NFOV 


NFOV 


prob 


N 50 = 0 . 75 


N 50 = 0 . 75 


N 50 = l . 50 


N 50 = 3 . 00 


N 50 = 6 . 00 


0.95 


9.74 km 


12.47 km 


10.87 km 


8.81 km 


6.19 km 


0.90 


10.21 


0.00 


11.24 


9.33 


6.78 


0.85 


10.48 


0.00 


11.49 


9.65 


7.18 


0.80 


10.70 


0.00 


11.69 


9.90 


7.50 


0.75 


10.89 


0.00 


11.86 


10.12 


7.79 


0.70 


11.06 


0.00 


12.00 


10.30 


8.02 


0.65 


11.20 


0.00 


12.13 


10.46 


8.25 


0.60 


11.34 


0.00 


12.26 


10.62 


8.47 


0.55 


11.48 


0.00 


12.39 


10.77 


8.68 


0.50 


11.62 


0.00 


12.52 


10.93 


8.87 


0.45 


11.76 


0.00 


0.00 


11.08 


9.08 


0.40 


11.90 


0.00 


0.00 


11.22 


9.30 


0.35 


12.05 


0.00 


0.00 


11.38 


9.51 


0.30 


12.20 


0.00 


0.00 


11.56 


9.73 


0.25 


12.38 


0.00 


0.00 


11.76 


9.98 


0.20 


12.59 


0.00 


0.00 


11.97 


10.27 


0.15 


0.00 


• 0.00 


0.00 


12.23 


10.57 


0.10 


0.00 


0.00 


0.00 


12.59 


11.01 


0.05 


0.00 


0.00 


0.00 


0.00 


11.65 



Figure 4.3 ACQUIRE typical output. Gunboat at Side aspect (range performance). N50 is a 
Johnson cycle criterion that corresponds to a 50% probability of target task (detection, 
recognition, classification, identification) 



D. WINEOTDA INPUTS 



WinEotda requires two kinds of input data: Meteorological and Site (Met) 



information, and Operations and Intelligence (Ops) data. The Met data includes the 



location, the surface weather forecast, and the thickness of the boundary layer. The data 



are entered using the Met form from the main screen (left comer). The Ops data includes 
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the electro-optical sensor_ID number, sensor height, viewing direction, targets, and 
backgrounds. The data are entered using the Ops form from the remaining three corners 
of the main screen. 

The following Figure shows the WinEOTDA Met data form. It will show up by 
using the input drop down menu button and selecting Meteorological Inputs. A detailed 
description of each parameter is enclosed in Appendix C. In order to differentiate the 




Figure 4.4 Met Data input form for WinEOTDA (Sample) 



two scenarios in this work, a set of meteorological data was used for each one based on 



the specific weather conditions of the previously referenced open field experiments at 



NAS Fallon, Nevada and off the coast of Monterey, California. 



41 





The available/required entries for the Met-input form can be selected from the 
three categories: Terminal Aerodrome Forecast (TAF), non-TAF code surface parameters, 
and radiosonde data. Since this is a laboratory work based on recorded data, the non-TAF 
surface parameters were used. 

The following Table shows the parameter set used for each one of the target 
scenarios. The Met data didn’t change when different sensors were implemented. 



Parameter 


Desert Environment 

Data from Reference 16 by 
Koch, Cynthia 


Maritime Environment 

Data from Reference 13 by 
Fu-Chau, Liu 


Surface Air Temperature 


32 F 


53 F 


Surface dewpoint Temperature 


28 F 


48 F 


Aerosol Index 


5 


12 


Battlefield Induced contaminants 


0 (none) 


0 (None) 


Visibility 


15 nm 


0 (None) 


Precipitation Index 


0 (None) 


0 (None) 


Rain Rate 


0 (None) 


0(None) 


Wind Speed 


15 knots 


8 knots 


Wind Direction 


300° 


298 ° 


Boundary layer 


400 hft 


400 hft 


Clouds 


Type: Low 
Height: 300 hft 
Amount: 8 


Type: Low 
Height: 300 hft 
Amount: 5 



Table 4.3 Meteorological data set used for each Environment 



The surface air temperature, surface dewpoint temperature, aerosol index, and 
visibility are used to compute an extinction coefficient. Extinction coefficient is used in 
ER, TV, and laser models. The extinction coefficients are used to compute the effect of 
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incoming radiation and the efficiency of signal transmission from the target scene to the 
sensor. An extended explanation for each parameter in this form is enclosed in Appendix 
C (pag. 108-110). 

The Ops data is entered using four different data entry forms: Sensor entry form, 
target entry form, target-location entry form, and background entry form. They are 
available by double clicking the descriptor picture in the main screen or from the drop 
down Input menu button. 

The following Figure 4.5 shows the target entry form used in this case: 





Target # 1 


Target #2 


Target 


|T72 T ank jJ 


(Frigate rj 


Heading 


fo“I jU jJ 


El -3 J 2d 


Position 






Op State 


j Exercised jJ 




Speed (kts) 


[10 jj J ,, ... 2d 





Target Elevation 
(ft MSL) 



w 



z I 




Figure 4.5 WinEOTDA Target Data Entry Form 
It shows a set of selected parameters for target #1: a T72, and also, the parameters 
for the target #2: a Frigate. Later, the target #2 was changed to a Gunboat for which the 
dimensions and characteristics are the same as the Research Vessel POINT SUR. A 
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detailed description of the parameters on this Figure 4.5 is enclosed in the Appendix C 
(pag. 107). 

The sensor entry form is represented with an aircraft on the main screen, double 
clicking on that picture will bring up the sensor form shown in Figure 4.6. 




Figure 4.6 WinEOTDA Sensor Data Entry Form 



Since this thesis deals with IR sensors the selection was made between 25 
different sensors. Also, a user-defined sensor # 150 was used with data obtained from the 
AN/AAA-36 physical parameters and theoretical values. The alphanumeric output shows 
the sensor data input for each run. See Appendix C (pag. 100-106). 
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The background data entry form is represented on the main screen by three 
rectangular boxes, with a generic picture of the selected background. Double clicking on 



each rectangular box brings up the following form: 



BACKGROUND DATA ENTRY FORM 



Background 81 



Water 


z J 


[Vegetation 


ISoil 


Snow 




Water 




Concrete 




Asphalt 




Swamp 




Rocky Field 





OK 



Figure 4.7 WinEOTDA Background Data Entry form 



The background affects the target signature in several ways. One is that the 
background temperature is a factor in the hot/cold spot selection. Another is the impact of 
background temperature on target mean temperature. So, the first-listed background is 
considered a primary background and plays a singular role. In computing the ground 
radiation that is reflected, TCM2 takes only the primary background into consideration 
[Ref. 1 l:p. A-6]. Since this fact could be a source of error, the author set the same target 
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T72, as target #1 and target #2 for the first run; and the Gunboat, as target #1 and target#2 
for the second run. 

The position of the target is defined by the latitude and longitude coordinates. 
The location of the target is used in the calculations of solar position, lunar position, and 
temperature. The form for this parameter is found at the bottom left comer, just below the 
target pictures, on the main screen. The two selected scenarios were very well located by 
their geographical coordinates, as shown in Table 3.3 in the previous chapter. 

When all parameters are set, the code WinEOTDA was run. So, when the 
execution is completed, an indication of a successful run is presented at the bottom of the 
main screen. Then, in the center of the main screen a box appears with the predicted 
maximum range for those specific parameters. 



E. WINEOTDA OUTPUTS 

There are three possible ways to get data output from WinEOTDA: 
Alphanumeric, graphic and tabular output. After the run is completed, the code 
automatically generates all three kinds of output. Also, the WinEOTDA results are 
presented in a grid on the main screen (pag. 96). In this case, just the maximum ranges 
for each target against the selected backgrounds are presented. The user has the 
possibility to change the units in the main screen grid-output by clicking on the label 
containing the units. The options are kilo-feet (kft), kilometers (km) and nautical miles 
(nm). 
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Tabular Output provides WinEOTDA ranges in an easy to read format. For the IR 
model, ranges are displayed for each target for both Narrow and Wide Fields of View. 
The range displayed is the longest detection range, factoring in all backgrounds and both 
MDT and MRT. 

The following Figure shows the table/tabular output for run # 8; the target is a 
Gunboat, the atmospheric conditions were set as a Maritime environment with strong 
continental influence and the user-defined sensor # 150 was used. 



IR EOTDA EXECUTION OUTPUT 



Time: 11/Q946Z Gunboat Gunboat 

Range (nm) Range (nm) 



View Direction 


NFOV 


WFOV 


NFOV 


WFOV 


' 0 00 


9.7 


6.6 


9.7 


6.6 


045 


12.2 


7.1 “ 


12.2 


7.1 


030 


12.6 


7.6 


126 


7.6 


_ '135 


12.4 


7.5 


12.4 


7.5 


M 180 


10.0 


6.5 


10 0 


6.5 




12.4 


7.5' 


12.4 


7.5 


270 


12.6 


7.6 


12.6 


7.6 


315 


12.2 


7.2 


12.2 


! 7.2 



0.0 — > No Value Computed 
-1.0 — > No Solution Possible 



Figure 4.8 WinEOTDA Tabular Output for run # 8 



The user can graph the contents of the table by using the right mouse button 
anywhere on the table. The author chose the 000 and 090 VIEW DIRECTION as the 
FRONT and SIDE aspect ratio in order to compare with the ACQUIRE results. 



47 



The Alphanumeric output is designed to permit creation of output products 
tailored to the user specific needs. All possible output products are created and are 
available for viewing. The code automatically produces output for the number of Fields of 
View (FOV) and the number of methods supported by the sensor. Since this output is too 
large a file to be shown here as a figure, it is enclosed as part of Appendix C (pag. 100- 
106). 

The Graphic output is designed to display a template with one to four plots. Each 
plot is a single graph defined by an x- and y-axis. The plot contents define which targets, 
backgrounds, and FOV are displayed. A sample of this kind of output is enclosed in 
Appendix C (pag. 1 14). 

F. RESULTS 

This chapter contains all procedures and intermediate results that were obtained 
through this work. This last section will show the final results and the comparison 
between these two codes. 

The Table 4.4 shows the predicted maximum ranges from ACQUIRE, assuming 
the required probability of detection based on the kind of task to be accomplished. Since 
WinEOTDA gives maximum range prediction only for detection, the remaining 
ACQUIRE results for classification, recognition and identification will be ignored for 
comparison purposes. However, a set of plots is presented as Figure 4.10 as a reference in 
order to show how ACQUIRE performs the different discrimination tasks. 
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As can be seen from the example in Figure 4.9, the results show some dependence 
on whether the sensor “sees” the target at Front or Side aspect. As was expected, the Side 
aspect presents greater values of range prediction, due of course to the increase in target- 



critical dimension. 



ACQUIRE: FRONT vs SIDE 




PROBABILITY OF DETECTION 



♦ NFOV(F) ■ NFOV(S) 



Figure 4.9 Relationship between Front and Side Aspect 



Figure 4.10 compares the ACQUIRE predicted range versus probability of 
detection for the T72 target (4.10 a) and the Gunboat target (4.10 b). In both cases, 
ACQUIRE produces larger values in NFOVd (detection mode). Also can be noticed that 
there is a difference for each performed task. 
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